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Abstract 

The European corn borer, Ostrinia nubilalis (Hiibner) (Lepidoptera: Crambidae), is a widely distributed and serious 
economic pest to corn production in the U.S. Genetic variability of 0. nubilalis was studied in 18 sub-populations in the 
upper Midwestern United States using amplified fragment length polymorphism. The relatively low Gsrvalues indicate 
that more variation exists within populations than between populations. High gene flow {Nm) values were indicated 
across the entire 0. nubilalis population; the lowest degree of gene flow was in the northern samples (Nm = 1.96) and the 
highest degree of gene flow was in the southern samples (Nm = 2.77). The differences observed in the respective regions 
(north vs. south) may be explained by the voltinism patterns (univoltine vs. multivoltine, respectively) of 0. nubilalis: 
southern multivoltine populations have opportunities for multiple matings for the duration of the year, further mix 
alleles. AMOVA results also indicated that most of the genetic variation was within sub-populations (~ 81% of total 
variation); less variation (~ 13%) was detected among populations within each of the three regions as designated for this 
study. However, the most striking and unexpected result was the low percentage of variation between all groups (=6%), 
further supporting implications of a high degree of gene flow. These results provide support for current requirements of 
refugia corn planting in Bt-corn management. These results also indicate that if resistance to Bt were to evolve in 0. 
nubilalis, quick action would be necessary to deter the rapid spread of the gene for resistance. 
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Introduction 

The European corn borer, Ostrinia nubilalis (Hiibner) 
(Lepidoptera: Crambidae), was introduced into the Un- 
ited States, and is a widely distributed pest that has 
proven to be a major biotic constraint for maize develop- 
ment and production. 0. nubilalis is known to be poly- 
phagous, attacking many herbaceous plants with stems 
large enough for the larvae to enter. Lewis (1975) repor- 
ted 223 plant species (both monocotyledon and dicotyle- 
don) on which the borers can develop. The existence of 
0. nubilalis in the United States was first reported by Vin- 
al (1917), however, 0. nubilalis is thought to have been in- 
troduced multiple times to North America in shipments 
of broom corn from Italy and Hungary into the eastern 
United States and Canada between 1909-1914 (Caffrey 
and Worthley 1927). Since its introduction, it has be- 
come one of the most destructive insect pests of maize in 
North America. 

0. nubilalis exhibits considerable genetic diversity. For ex- 
ample, voltinism differences between populations in 
North America were recognized shortly after the insect 
was discovered. Voltinism associated with diapause is an 
inherited characteristic modified by environmental 
factors. The wide geographic distribution exposes 0. nu- 
bilalis to ecological conditions that differ in photoperiod, 
temperature, host plant availability, and growing season 
length (Calvin et al. 1991). Populations can be character- 
ized as univoltine, bivoltine, and multivoltine (Showers et 
al. 1989). Ikten (2002) found voltinism to be genetic in 
origin, sex-linked, and controlled by a few loci. Further, 
Ikten (2002) found voltinism displayed a response to short 
photoperiods whereby 0. nubilalis could adapt quickly to 
local conditions. As a result, historically bivoltine popula- 
tions can become univoltine by a simple drop in temper- 
ature during critical days of diapause giving the species 
the flexibility to take advantage of the full growing season 
depending on altitude and latitude. 

While pesticides have helped control insect pests of food 
and fiber, pesticide resistance will occur for any pest con- 
trol product with high selection pressure that removes 
susceptible individuals from the population and allows 
only those individuals that possess resistance genes to re- 
produce with no other intervention (Clark and Yamagu- 
chi 2002). Biotechnological developments, such as the 
isolation of genes with insecticidal properties and the 
ability to insert those genes into sexually incompatible 
species, have allowed scientists and producers to diminish 
insect damage to plants. An effective application of bio- 
technology has been the expression of bacterial insect- 
icidal toxins from Bacillus thuringiensis Berliner (Bt) in sev- 
eral plant species (i.e. maize, cotton) where they can 
provide control of lepidopteran insect pests such as 0. nu- 
bilalis. With current demands by producers for Bt crops, 
there is potential for the development of resistance to 
plant mediated Bt toxins. It is important to develop 
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better understanding of the insect's genetic structure, ge- 
netic variation, and gene flow that can provide the basis 
for improvement and changes in current management 
strategies for insect control and resistance management. 

A great deal of research has been devoted to the biology 
and behavior of 0. nubilalis, driven by the need to de- 
termine the affects of 0. nubilalis injury to crops, particu- 
larly corn, %ea mays, to assess economic damage and to 
develop management strategies (Hunt et al. 2001). 
However, limited information is available about 0. nu- 
bilalis genetic variation, gene diversity, and gene flow. 
Previous studies comparing allozymes (Harrrison and 
Vawter 1977; Carde et al. 1978; Cianchi et al, 1980; and 
Glover et al. 1990) of pheromone and voltinism races 
found limited genetic variation while Marcon et al. 
(1999) using PCR-RFLP found similarity among 0. nu- 
bilalis populations. From a resistance management stand- 
point, this has major implications for selection pressure in 
localized areas resulting in a higher likelihood of resist- 
ance to pesticides including plant-mediated Bt. However, 
Pornkulwat et al. (1998) and Saldanha (2000) found evid- 
ence of genetic variability among 0. nubilalis populations 
using random amplified polymorphic DNA (RAPD) 
markers. Bourguet et al. (2000) also studied gene flow of 
French populations of 0. nubilalis and found a high and 
homogenous gene flow. More recently, Coates et al. 
(2004) found significant genetic differentiation between 
Atlantic coast and Midwestern United States samples. 
But Bazin et al. (2006) indicated that there does not ap- 
pear to be a detectable correlation between mtDNA 
polymorphism and species abundance; mtDNA may re- 
flect the time of the last selective sweep rather than the 
demographics or population history of the organism and, 
concluded that mtDNA diversity is unpredictable making 
it not useful for biodiversity studies. 

Genetic variation, governed by natural selection through 
the interaction of genetic forces and changing environ- 
ments through space and time, provides the basis for 
evolutionary change. Genetic differentiation between 
populations largely depends on the interacting balance 
between gene flow, genetic drift, and natural selection 
(Futuyma and Peterson 1985). Although natural selection 
acts direcdy on phenotype, it is a major factor causing ge- 
netic differentiation at the protein and DNA level. Stud- 
ies of genetic variability used to infer gene flow could 
help determine if the entire 0. nubilalis population is truly 
a single interbreeding population in which intense gene 
flow can occur. 

The AFLP technique has proven to be a reliable tool to 
generate highly polymorphic molecular markers used to 
study genetic divergence in insect populations such as: 
Spodoptera frugiperda (McMichael and Prowell, 1999; Suin- 
aga et al. 2004; Clark et al. 2006), Lymantria dispar 
(Reineke et al. 1999), and Bemisia tabaci (Cervera et al. 
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2000). The AFLP technique combines the reliability of 
RFLPs and power of PGR in a single technique. Garcia 
et al. (2004) compared RAPD, RFLP, AFLP, and SSR 
markers and concluded that AFLP markers were the best 
choice for evaluating the diversity and genetic relation- 
ships between tropical maize lines without requiring pre- 
vious knowledge of any DNA sequencing. Alamalakala 
(2002) found that AFLP consistently distinguished 
between European and North American populations of 
0. nubilalis providinga good starting point for further in- 
sect genetic studies using AFLPs for 0. nubilalis. 

The objective of this study was to measure genetic vari- 
ation within and between sub-populations to infer genetic 
diversity and gene flow for 0. nubilalis. Current informa- 
tion on gene flow is limited, but knowledge of the degree 
of genetic variation and gene flow of 0. nubilalis is imper- 
ative to develop effective management strategies for crop 
protection. 



Materials and Methods 

Sample collection 

Twenty specimens were collected from 17 sub-popula- 
tions (locations) in eight states in the upper Midwest and 
one location from Louisiana (LA) (Table 1; Figure 1). 
Specimens were frozen in an ultra-low freezer to kill 
them quickly, maintain DNA integrity, and for long term 
storage. All samples, except from Louisiana, were over- 
wintering 5 instars collected in the fall from maize 
which allowed all ecotypes (univoltine, bivoltine, mul- 
tivoltine) of 0. nubilalis to be present. The LA adult 
samples were collected from rice in mid summer. 

DNA isolation and quantification 

DNA was isolated from 10 of the 20 individually frozen 
insect samples using a modified Black and Duteau (1997) 
CTAB extraction protocol. Prior to homogenization, the 
gut was removed from the larvae; and the head, wings, 
and abdomen were removed from adult insects from Mer 
Rouge, LA. The procedure was modified with the addi- 
tion of a chloroform/ phenol extraction step prior to eth- 
anol precipitation. After removal of ethanol, the pelleted 



Table I. Code, collection site, and year of collection of 0. nubilalis populations. 



Code 


Collection Site 


State 


Year 
Collected 


Elevation (m) Above Sea 
Level 




Northern Sample 
Sites 








CL 


Clay County 


SD 


2003 


372 


SS 


South Shore 


SD 


2003 


568 


EL 


Ellendale 


ND 


2003 


444 


FR 


Fargo 


ND 


2003 


326 


WM 


Willmar 


MN 


2003 


329 


CO 


Concord 


NE 


2003 


349 




Western Sample 
Sites 








NP 


North Platte 


NE 


2003 


878 


OK 


Oshkosh 


NE 


2003 


1073 


GR 


Gordon 


NE 


2003 


1083 


WR 


Wray 


CO 


2003 


1 134 


TR 


Torrington 


WY 


2003 


1257 


AL 


Alliance 


NE 


2003 


1207 




Southern Sample 
Sites 








MD 


Mead 


NE 


2003 


367 


CC 


Clay Center 


NE 


2003 


546 


AB 


Auburn 


NE 


2003 


314 


TK 


Topeka 


KS 


2004 


329 


CM 


Cameron 


MO 


2004 


316 


LA 


Mer Rouge 


LA 


2003 


29 
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DNA was dissolved in 50 ^1 IX TE buffer (lOmM Tris- 
HGLpH 8.0; 0.1 mM EDTA). 




Figure I. Sample distribution of 0. nubilalis sub-populations 
across the upper Midwestern United States and in Louisiana. 



Quantification to a known concentration of X DNA (Life 
Technologies, www.invitrogen.com) was performed by 
loading 1.0 fA of DNA mixed with 1.0 fA tracking dye to 
a 1% agarose gel that was electrophoresed for 30 min. 
Ethidium bromide stained gels were then visualized using 
Advanced Quantifier version 4.0 

(www.genomicsolutions.com). DNA samples were diluted 
with TE buffer if concentrations were >22.5 ng/ fA. 

AFLP-PCR 

The AFLP procedure was achieved using 3 steps: 1) 
DNA template preparation; 2) DNA template preampli- 
fication and; 3) AFLP Selective amplification. Template 
preparation and AFLP assays were performed using a 
modified protocol of Vos et al. (1995). Infrared labeled 
(IRD-700) (LI-COR, Lincoln, NE, USA) EcoRl primers 
were used in the polymerase chain reaction. 

Template preparation 

Approximately 160 ng of genomic DNA in 7 fA of 
volume was incubated with restriction endonucleases 
_£coRI and Msel (New England Biolabs, www.neb.com) 
for 2.5 hrs at 37°C in a total volume of 12.5 [A which 
contained 1.25 fA of 10X One-Phor-All buffer 
(Amersham Pharmacia biotech, www.apbiotech.com), 



0.125 fA of 10 U//A Msel enzyme (1 .25U/ reaction), 
0.0625 fA of 20 U/fA EcoRl enzyme (1.25 U/reaction), 
0.125 fA of 10 mg/ml BSA (bovine serum album) (New 
England Biolabs), and autoclaved nanopure water to 
make up the volume to 12.5 fA. The restriction fragments 
were then incubated for 1 1.5 hrs at ~ 25°C with 5 fA lig- 
ation mixture containing 0.15 fA of T4 DNA ligase, 0.5 
fA of 10X T4 ligase buffer (New England Biolabs), 0.5 fA 
of 5 pmoles/^1 EcoRl adapter, 0.5 fA of 5 pmoles///l 
Msel adapter (Operon Technologies, www.operon.com), 
and 3.35 fA autoclaved nanopure water. After ligation, 
the template was diluted 8-fold. Double stranded ad- 
apters for ligation to the corresponding cut made by the 
restriction enzyme were prepared by incubating equimol- 
ar amounts of both adapter strands for 10 min at 65°C, 
10 min at 37°C, 10 min at 25°C and stored at -20°C. 

Preamplification of DNA template 

PCR amplification of the ligated material consisted of 20 
cycles (30s at 94°C, 1 min at 56 C, and 1 min 72°C) us- 
ing 1.0 fA of diluted template from adapted DNA diluted 
8-fold with IX TE buffer (lOmM Tris-Cl, 0.1 mM EDTA 
[pH8.0]), 8 fA pre-amp primer mix II (contained two oli- 
gonucleotide primers, one corresponding to the EcoRl 
adapted ends and one corresponding to the Msel ends; 
Life technologies, www.operon.com), 1.0 fA of 10X PCR 
buffer containing 15mM MgCl2 and 0.25 /A of 5 U/fA 
AmpliTa<7 DNA polymerase (1.25 U/reaction) (Applied 
Biosystems, www.appliedbiosystems.com). The oligonuc- 
leotide primers in the pre-amp primer mix II are comple- 
mentary to the adapter/ restriction site with Msel primer 
containing one selective nucleotide M(N+1) primer and 
iicoRI primer containing no selective nucleotide E(N+0) 
primer (Table 2). Products from the preamplification re- 
actions were diluted 20-fold with autoclaved nanopure 
water and were used as the template for selective ampli- 
fication reactions. 

Selective PCR amplifications 

After screening 24 primer combinations, five primer 
combinations (Table 3) were selected for AFLP analysis. 
Reaction volumes of 10.5 fA containing 2.0 fA of diluted 
preamplified template, ~ 5.0 fA of autoclaved nanopure 
water, 1.2 fA of 10X PCR buffer containing 15 mM 
MgCl2, and 0.06 fA of 5 U/fA AmpliTaq polymerase 
(Applied Biosystems), 2.0 fA of Msel primer (M+l+2) 
(6.7 ng/fA dNTPs) (Life Technologies), and 0.2-0.5 fA of 
1.0 pinoles//*! IRD-700 labeled EcoRl primer (LI-COR, 
www.licor.com). Selective amplifications were conducted 
in a DNA thermal cycler 9600 (Applied Biosystems) using 
the following Touchdown PCR program: 1 cycle of 30s 
at 94° C, 30s at 65° C and 1 min at 72° C; 12 cycles of 
30s at 94° C, 30s at 65°C, (subsequently lowering the an- 
nealing temperature by 0.7°C per cycle) and 1 min at 
72° C; 23 cycles of 30s at 94° C, 30s at 56° C, 1 min at 
72° C. After selective amplifications, reactions were 
stopped by adding 2.5 fA blue stop solution (LI-COR). 
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Table 2. Oligonucleotide adapters and primers used for AFLP analysis. 



Primer ID 


Primer Type 


Sequence (5 -3 ) 


EcoRI-F 


Forward Adapter 


CTCGTAGACTGCGTACC 


EcoRI-R 


Reverse Adapter 


AATTGGTACGCAGTCTAC 


Msel-F 


Forward Adapter 


GACGATGAGTCCTGAG 


Msel-R 


Reverse Adapter 


TACTCAGGACTCAT 


E(N+0) 


EcoRI Preamp Primer 


GAC TGC GTA CCA ATT C 


M(N+I) 


Msel Preamp Primer 


GAT GAG TCC TGA GTA AC 


E(N+0+3) 


EcoRI Selective Primer 


GAC TGC GTA CCA ATT C NNN* 


M(N+l+2) 


Alsel Selective Primer 


GAT GAG TCC TGA GTA A CNN* 



N represents the position of (+3 and +2) selective nucleotides in the primer and may be either A,T,G, 
or C. 



Table 3. Selective Primers used for AFLP analysis their associated amplified loci and range of fragment 
sizes for O. nubilalis. 



Primer Pair + Selective 
Ext. 


Number of Polymorphic 
Loci 


Fragment Size Range 
(bp) 


EcoRI + ACA/Msel + CTC 


33 


57-21 1 


EcoRI + ACT/Msel + CTG 


27 


71-242 


EcoRI + ACT/Msel + CAT 


26 


72-239 


EcoRI + ACT/Msel + CAA 


36 


70-34 1 


EcoRI + ACA/Msel + CTT 


42 


66-262 


Total 


164 


57-341 


Average/Primer Pair 


33 





Samples were then denatured at 95 C for 3 min and 
flash cooled on ice immediately before loading the poly- 
acrylamide gel for electrophoresis. One fi\ of the sample 
was electrophoresed through 

KB plus 

6.5% ready-to-use 

gel matrix (LI-COR), and the infrared fluorescent bands 
were then detected by a laser scanning system LI-COR 
Model 4200S-2. 

Scoring AFLP data 

Sizes of AFLP fragments were estimated using an 
IRD-700 labeled 50-700bp ladder. Fragments were 
scored using SAGA MX version 3.2 (LI-COR); a soft- 
ware program designed specifically for scoring AFLP 
data. The data were then converted to numerical data (Is 
for band presence, 0s for band absence) that allowed the 
identification of polymorphic markers. There were five 
gels for each primer pair containing four populations 
randomly picked for a total of 25 PAGE gels. All gels 
were viewed for clarity prior to scoring to establish scor- 
able bands at a given fragment size for all populations. 
All bands consistently recognized by SAGA Generation 2 
Software Version 3.2 (LI-COR), both monomorphic and 
polymorphic, were scored and included in the binary 
data matrix. 



Coefficient of variation analysis 

Analysis was conducted using DBOOT version 1.1 
(Coelho 2001) to assess the appropriate number of poly- 
morphic loci required for acceptable precision for genetic 
analyses. Bootstrapping analysis using the simple match- 
ing coefficient (Kosman and Leonard 2005) was conduc- 
ted with 1000 permutations to assess the strength of the 
molecular markers. A comparison was then made 
between the coefficient of variation values against the 
number of markers scored and was then plotted. 

Genetic diversity and gene flow of O. nubilalis 

Data were analyzed using POPGENE version 1.32 (Yeh 
and Boyle 1997) using a dominant marker data set (164 
markers) assuming Hardy-Weinberg equilibrium. 
Samples were grouped into three geographical regions 
due to proximity. Analyses were observed at three levels 
(1) individuals; (2) regions; (3) and whole populations. 
The percent (%) polymorphism, genetic diversity or het- 
erozygosity (H), Gst, an d gene flow estimation (Not) 
were then assessed within and between all populations. 
Individual populations were analyzed for genetic di- 
versity (H) for each sub-population as per Nei (1973). 
Gst values generated by POPGENE were given for the 
whole population and were expressed as the gene 
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diversity of the single population subtracted from the 
gene diversity of the total population divided by the gene 
diversity of the total population (GsT — Ht ~~ Hs / Ht). 
Gene flow was estimated from GsT values and is ex- 
pressed as (Nm) = 0.5(1 - Gst)/Gst (McDermott and 
McDonald 1993). 

Mantel tests comparing genetic dissimilarity 
and geographical distance 

A permutation test (Mantel 1967) was conducted to test 
the relationship between geographical distance and the 
measure of genetic dissimilarity (1 -similarity), to further 
verify the estimated levels of gene flow found in the 
POPGENE analysis. The Mantel permutation tests were 
conducted using the MxCOMP module in the NTSYSpc 
ver. 2.1 (Rohlf 2002) software program with 1000 per- 
mutations. Tests were first conducted on the full data set 
and then separately with western samples, northern 
samples, and the southern samples as data sets. 

Analysis of molecular variance 

Analysis of molecular variance (AMO VA; Excoffier et al. 
1992) was conducted with ARLEQUIN 2.0 (Schneider et 
al. 2000). In this analysis, total variance of the AFLP data 
set was partitioned at three hierarchical levels: (1) an 
among-population component; (2) a regional or six sub- 
population component; (3) and a within-population com- 
ponent. Unlike the calculations used for Nei's GsT val- 
ues, the AMOVA partitions the variation according to 
correlations among genotypes rather than gene frequen- 
cies due to the dominant expression of AFLP markers. 
The significance of the three variance components was 
tested using 1000 random permutations. A two-part 
AMOVA analysis was conducted to test genetic diver- 
gence (Fst) a s a factor of variation among individuals 
within a given population and between populations. 

Principal component analysis (PCOA) 

Principal component analysis was used to test genetic 
isolation between sampled populations using NTSYSpc 
ver. 2.1 (Rohlf 2002). This two-dimensional analysis, an 
Eigen plot, helps identify patterns in the data and enables 
one to graph the results highlighting the similarities and 
differences. 

UPGMA cluster analysis 

An Unweighted Pair Group Method with Arithmetic 
Mean (UPGMA) consensus cluster analysis as outlined by 
Sneath & Skoal (1973) was conducted using NTSYSpc 
ver. 2.1 (Rohlf 2002) on all 18 sub-populations of 0. nu- 
bilalis to illustrate genetic similarity. Bootstrap analysis, 
using BOOD-P software ver. 3.1 (Coelho 2001), was 
used to test the reliability of the dataset with 10,000 
pseudoreplicates. 



Results 

Number and size of AFLP loci observed 

The five primer combinations used with DNA from 18 
sub-populations and 1 80 individual samples of 0. nubilalis 
produced 164 amplified loci, averaged 33 loci per primer 
combination, and fragment size ranged from 57—341 bp 
(Table 3). Polymorphism within the 18 sub-populations 
was high, averaging 83% (136 out of 164 were poly- 
morphic) and ranging from 72% to 94% (Table 4), indic- 
ating that AFLP provides polymorphic markers for study- 
ing genetic variation. 



Coefficient of variation of AFLP loci 
(markers) explored 

The correlation of the coefficient of variation and the 
number of molecular markers examined for the sampled 
population is important for determining robustness in ge- 
netic variability studies because either a large number of 
loci or a large number of individuals should be examined 
(Hoelzel 1995). Analysis using the DBOOT program 
showed that the high number of markers used decreased 
the coefficient of variation to the point that all but 6.7% 
of the variation in the population was explained (Figure 



Table 4. Heterozygosity and polymorphism (%) for 18 sub- 
populations of O. nubilalis. Individual heterozygosity (HS) indic- 
ates a lack of homogeneity. 



Location 


State 


Bands 


(%) 


Gene Diversity 
(HS) 


Mead 


NE 


1 18 


72.0 


0.2370 


Clay Center 


NE 


129 


78.7 


0.2879 


Concord 


NE 


139 


84.8 


0.3054 


North Platte 


NE 


154 


94.0 


0.3757 


Oshkosh 


NE 


129 


78.7 


0.3083 


Alliance 


NE 


137 


83.5 


0.2956 


Gordon 


NE 


137 


83.5 


0.3307 


Auburn 


NE 


144 


87.8 


0.31 17 


Clay County 


SD 


142 


86.6 


0.3332 


South Shore 


SD 


143 


87.2 


0.3 1 68 


Ellendale 


ND 


1 18 


72.0 


0.2651 


Fargo 


ND 


135 


82.3 


0.3039 


Willmar 


MN 


141 


86.0 


0.3140 


Wray 


CO 


123 


75.0 


0.2681 


Torrington 


WY 


135 


82.3 


0.3222 


Topeka 


KS 


144 


87.8 


0.3209 


Cameron 


MO 


153 


93.3 


0.3705 


Mer Rouge 


LA 


133 


81.1 


0.3029 


Average 




136 


83.1 




Total 
Population 




164 


100.0 
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Figure 2. The coefficient of variation of 0. nubilalis populations using AFLP molecular markers with 1 000 bootstrap replications. 
Greater than 90% of the variation in the 0. nubilalis population is explained by these AFLP markers. 



2). Therefore, 164 markers were sufficient for unbiased 
genetic analyses. 

Genetic diversity, Gst values, and gene flow 

Genetic diversity, or population heterozygosity (H), val- 
ues calculated were high from all sub-populations, ran- 
ging from 0.2370 at Mead, (NE) to 0.3757 at North 
Platte, (NE); overall gene diversity averaged 0.3094 
(Table 4). Analysis across all populations revealed a high 
genetic diversity value of 0.4121 (Table 5). The Gst 
value among all sampled populations was low, 0.1652, 
indicating a high degree of within population variation 
(83%) and low variation among populations (16.5%). Not 
values were high in all three regions (Table 5), suggesting 
that events, such as migration, provide opportunity for 
interbreeding and gene flow. 

Permutation tests comparing genetic dissim- 
ilarity and geographical distance 

Results indicated no evidence of correlation between ge- 
netic distance and geographical distance (data not 
shown). There was no correlation for the western samples 
(r = -0.3765, t = 1.484, P = 0.9312), northern samples (r 
= 0.2937, t = 1.1087, P = 0.8662), southern samples (r = 
0.2546, t = 0.9671, P = 0.8332), nor for the full data set 
combined (r = -0.0146, t = -0.1223, P = 0.4513). These 
results further imply there is not genetic isolation in this 



portion of 0. nubilalis range and agree with results of 
Coates et al. (2004). 

Analysis of molecular variance (AMOVA) 

The AMOVA analysis showed that approximately 81% 
of the variation in the data set was from genotypic vari- 
ation within populations (Table 6). Only 13% of the vari- 
ation could be attributed to differences among popula- 
tions within regions while the remaining 6% was due to 
the variation among regions. The genetic divergence 
(Fst) value for 18 sub-populations was 0.1907 indicating 
low differentiation between populations. These results 
were concordant to those obtained using POPGENE 
(Table 5). Similar results were reported by Coates et al. 
(2004); this low level of population differentiation sup- 
ports the possibility of gene flow across this portion of the 
0. nubilalis range. 

Principal component analysis (PCOA) 

Principal component analysis for these 18 sub-popula- 
tions of 0. nubilalis revealed four distinct clusters (Figure 
3). Cluster 1 consisted of mixed populations from several 
states that may represent areas of transition from univolt- 
ine to bivoltine ecotypes, cluster 2 represented univoltine 
samples from South Dakota and North Dakota, cluster 3 
represented samples typical of the univoltine ecotype, 
and cluster 4 represented samples that were less similar 
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(more scattered on the graph) that are southern popula- 
tions probably representing the transition from bivoltine 
to multivoltine ecotype. 



UPGMA cluster analysis 

Consensus averages yielded a dendrogram with similarity 
coefficients ranging from 65% to 74% for all 18 sub-pop- 
ulations of 0. nubilalis, indicating a high degree of 



variability within locations (Figure 4). Bootstrap values 
ranged from 26 — 100% but a majority of nodes had 
bootstrap values exceeding 70%. Most sub-populations 
clustered with others in close geographic proximity or 
with similar climate. For example, samples from the 
western region (higher altitude) and from the northern 
region (higher latitude), with a similar level of growing 
degree days and suspected to consist primarily of a uni- 
voltine ecotype, were more genetically similar than 
samples from southern locations. Those nodes showing 
lower bootstrap values (AB, CC, CM and CO) indicate 
weak association and could represent areas of transition 
with a mixture of ecotypes (particularly bivoltine and 
multivoltine) as is also indicated with the PCOA analysis 
(Figure 3). As reported by Pornkulwat et al. (1998), uni- 
voltine ecotypes, using current techniques, may be more 
easily identified than bivoltine or multivoltine ecotypes, 
indicating that gene flow may be somewhat restricted for 
the univoltine ecotype. 

Discussion 

Molecular genetic tools have proven successful at detect- 
ing genetic variability. Pornkulwat et al. (1998) used 
RAPD markers that were able to distinguish multivoltine 
from univoltine and bivoltine ecotypes. Saldanha (2000) 
used PvAPD-PCR to distinguish between local popula- 
tions of 0. nubilalisand found a large genetic group con- 
sisting of univoltine, bivoltine, and multivoltine ecotypes 
in Nebraska. However, disadvantages of RAPDs are ap- 
parent and therefore results can be arguable. 
Alamalakala (2002) showed that using AFLPs allows the 
analysis of large numbers of loci while repeatability is 
greatly increased over RAPDs. Our AFLP studies evalu- 
ated numerous loci showing a very high degree of poly- 
morphism between the 1 8 sub-populations of 0. nubilalis 
studied. 

Understanding population structure of 0. nubilalis will 
provide critical base-line information for developing sus- 
tainable management strategies. Populations may either 
exist as a single panmictic unit, as if in Hardy-Weinberg 
equilibrium, where all individuals have the equal oppor- 
tunity to randomly mate, they may occur as a series of 



Table 5. Pair wise comparisons of Nei's coefficient of gene 
differentiation (Gsr) between three regions of O. nubilalis pop- 
ulations and estimates of gene flow (Nm). 



Location 


Heterozygosity 


Gst 


Nm 


Northern Sample 
Sites 


0.3004 


0.203 1 


1.9617 


Clay County 


SD 








South Shore 


SD 








Ellendale 


ND 








Fargo 


ND 








Willmar 


MN 








Concord 


NE 








Western Sample 
Sites 


0.3393 


0.1719 


2.4081 


North Platte 


NE 








Oshkosh 


NE 








Gordon 


NE 








Wray 


CO 








Torrington 


WY 








Alliance 


NE 








Southern Sample 
Sites 


0.3809 


0. 1 529 


2.7696 


Mead 


NE 








Clay Center 


NE 








Auburn 


NE 








Topeka 


KS 








Cameron 


MO 








Mer Rouge 


LA 








All Populations 


0.4121 


0.1652 


2.5261 



Table 6. Hierarchical analysis of molecular variance (AMOVA) for 18 sub-populations of O. nubilalis. The 
majority of the variation is within a given population. 



Source of Variation 


d.f. 


Sum of 
Squares 


Vaiance 
Components 


Percentage of 
Variation 


Among Groups 


2.0 


386.3 


2.0 


6.0 


Among Populations Within 
Groups 


15.0 


1 070.4 


4.4 


13.1 


Within Populations 


162.0 


4424.6 


27.3 


80.9 


Total 


179.0 


5881.3 


33.8 
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Figure 3. Principal Component Analysis (PCOA) of 1 8 sub-populations of 0. nubilalis. Four distinct clusters were evident pos- 
sibly indicating some genetic isolation from other clusters. 



small sub-populations which are isolated from one anoth- 
er (island model), or they may exist as a continuous popu- 
lation where genetic exchange takes place by geographic- 
ally proximate individuals and there is some genetic isol- 
ation by distance (Slaktin and Barton, 1989). Frankham 
et al. (2002) indicated that a Gst value of 0.15 indicates 
differentiation among fragments. At the population level, 
our study indicated a Gst value of 0.1652. Therefore, 0. 
nubilalis in the Midwest may not consist of a continuous 
population as some differentiation may be taking place, 
as discussed by Coates et al. (2004), but this is far re- 
moved from representing an island model. Studies by 
Bourguet et al. (2000) also concluded there is a high de- 
gree of gene flow and no differentiation between north- 
ern (univoltine) and southern (bivoltine) populations in 
France: observed and expected heterozygosities were al- 
most identical 0.19 to 0.33 and 0.23 to 0.32 respectively 
based on six polymorphic allozyme markers. 

Like the Gst analysis, AMOVA results indicated that a 
majority of the variation detected was within sub-popula- 
tions. Only a small amount of variation was due to vari- 
ation among the three regions studied (north, south, and 
west). However, the most striking and unexpected result 
was the low percentage of variation among groups, 



further enhancing the implication of a high degree of 
gene flow possibly through migration events. The degree 
of variation within or among populations was variable. 
For many species, variability may be greater between 
than within populations (Roderick 1996; Clark et al. 
2006). Local selection pressures and barriers could 
provide for genetic isolation resulting in a high degree of 
genetic drift. However, previous studies showed a higher 
degree of variation within populations than between pop- 
ulations for many insects (Coates and Hellmich 2003, 
Juan et al. 2004, Timmermans et al. 2005). Results from 
our study examining 0. nubilalis populations from the 
Midwestern U.S. using AFLPs show a larger degree of 
genetic variation than previously reported but is similar 
to results found by Bourguet et al. (2000) with 0. nubilalis 
in Europe. 

The numerical differences observed in Not between re- 
spective regions (north vs. south) and respective UPGMA 
clustering may be explained by the voltinism patterns 
that prevail in those regions. The opportunity for mul- 
tiple mating during the year in southern regions could 
result in high genetic variability. Another possible ex- 
planation of lower indications of gene flow in the north- 
ern region could be that 0. nubilalis living on alternate 
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Figure 4. UPGMA Dendrogram of 18 sub-populations of 0. nubilalis. Bootstrap values are at each node. 



hosts, such as Solanum and Amaranthus may be unable to 
mate with individuals from maize due to different emer- 
gence patterns of 0. nubilalis from the alternate hosts 
compared to those from maize. Or, as suggested by 
Coats et al. (2004), there may be genetic differences 
between sympatric uni- and bivoltine ecotypes. Coats et 
al. (2004) evaluated individuals using a portion of 
mtDNA whereas we examined 164 loci from random 
locations across the genome of 180 individuals. It would 
be beneficial to conduct similar studies over a broader 
range using a greater number of individuals for genetic 
analyses using AFLP that provides the polymorphic 
markers necessary to better address the questions of ge- 
netic differentiation in 0. nubilalis. 

Showers et al. (2001) conducted mark and recapture 
studies to track movement of 0. nubilalis after release. 
They found 37-52% of 0. nubilalis males recaptured flew 
800 m or more, 8—1 1% flew 3.2 km or more, while some 
collections of marked males and females occurred at dis- 
tances greater than 40 km. However, a very large 



percentage of the 0. nubilalis released were never recap- 
tured (over 99.6%). It is possible that some of these indi- 
viduals traveled much further than placement of the re- 
capture pheromone traps. Showers et al. (1995) reported 
similar findings while studying potential movement on 
surface airflow of 0. nubilalis. Although indirect measures 
of gene flow are tenuous (Whitlock and McCauley 1 998, 
Hey and Machado 2003) our results, along with previous 
results (Showers et al. 1995, 2001) illustrating that 0. nu- 
bilalis were capable of dispersing at least 32 km/yr, would 
support the concept of gene flow within a continuous 
population. Results from these ecological and genetic 
studies provide support for current refuge requirements 
relative to the use and management of Bt-corn. Also, 
these results indicate that if insect resistance was to evolve 
in 0. nubilalis, quick action will need to be taken to defer 
the spread over large geographic areas. 
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